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ABSTRACT 


In  certain  applications  involving  reflection  from  smooth  convex  surfaces, 
the  conventional  geometrical  optics  method  (GO)  is  not  sufficiently  accurate. 

In  this  paper  a  generalization  of  the  conventional  GO  formulas  is  presented  and 
closed-form  expressions  are  derived.  The  derivation  of  the  generalized  geo¬ 
metrical  optics(GGO)  formulas  is  based  on  the  Fresnel  transformation  between 
two  planes,  one  of  which  contains  the  spectral  point,  while  the  other  passes 
through  the  observation  point.  The  method  developed  herein  improves  the 
accuracy  of  GO  without  paving  the  high  numerical  cost  of  the  physical  optics 
integration.  It  is  important,  in  particular,  for  computing  the  reflection 
characteristics  of  numerically  specified  surfaces. 
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High-frequency  scattering  of  an  electromagnetic  field  from  a 
smooth  convex  conducting  surface . 

Axial  and  paraxial  rays  in  a  doubly  divergent  pencil . 

Coordinate  system  pertinent  to  the  continuation  of  a  field 
distribution  from  one  plane  to  another . 

Aperture  planes  3q(z  •  Zq)  and  o^(z  •  z^)  normal  to  the  reflected 

rav  at  the  specular  point  A  and  the  observation  point  P, 
respectively . 


The  transverse  axes  (x,v)  make  an  angle-”.-  with  the  principal 
axes  (x_,y_) . 


1.  INTRODUCTION 


Most  of  che  dif f ractlon-theoretic  approaches  co  high-frequency 
scattering  from  a  smooth  conducting  reflector  involve  integration  over 
the  current  induced  on  the  reflector  by  a  given  incident  field  [1].  The 
most  straightforward  of  these  approaches  is  the  Physical  Optics  (PO) 
approach,  wherein  the  currents  on  che  reflector  are  approximated  by  the 
currents  calculated  from  the  theory  of  Ceometrical  Optics  (GO).  If  no 
characteristics  of  che  integrand  C3n  be  found  that  will  permit  simplifica¬ 
tion  of  the  PO  integral,  it  is  often  necessary  to  employ  a  "brute  force" 
integration  which  turns  out  to  be  lengthy  and  laborious.  However,  if  the 

reflector  surface  is  doubly  convex,  i.e.,  non-focusing,  the  PO  integral 

can  be  as-mpcoticallv  evaluated  at  high  frequency  using  the  method  of 
scationarv  phase.  The  asvmptotlc  physical  optics  yields  closed-form 
expressions  that  are  identical  to  those  given  by  the  ray  method. 

However,  in  certain  applications,  the  ray  method  result  is  not 
sufficiently  accurate,  and  the  higher  order  corrections  to  GO  are  necessary. 
The  classical  wav  to  obtain  higher  order  GO  terms  is  to  expand  the  electro¬ 
magnetic  field  in  a  high-frequency  asvmptotlc  series  of  the  form 

E  ■  exp  (-Jks(r)l  (jk)  *  e  (r)  (1.1) 

w  IB 

B“0 


where  s  is  the  phase  function  and  {e  }  are  amclitude  vectors.  Using  ray 
techniques,  Lee  (2,3)  was  able  to  obtain  explicit  formulas  for  calculating 
the  first  two  orders  in  (1.1)  of  the  field  reflected  from  a  conducting 
surface  in  terms  of  the  incident  field  and  the  surface  properties  at  the 


specular  point  A.  (See  Fig.  1.)  However,  continuation  of  the  amplitude 
vectors  e  from  the  specular  point  A  (z»z.)  to  the  observation  point  P 
entails  the  use  of  recursion  relations  of  the  form: 

e0U>  -  *0U0>  DF(z , zQ)  (1.2a) 

e  (z)  •  e  (z  )  DF(z,z  ) 

a  a  0  0 

z 

-  1/2  /  DF(z.z’)  7‘e  ,(z')  dz'  m  -  1,  2,  3.  ...  (1.2b) 

*  !D-  L 

2o 


where  DFiz.z^)  in  the  divergence  factor  (to  be  explicitly  defined  later) 

and  a  (z)  and  7‘e  iz)  are  given  bv 
m  *a 


e  (z)  - 

7‘e  (z) 

a 


(x,y,z) ] 

x»y«0 


-  (7~^(x,y.z)  ] 

x"v»0 


(1.3a) 

(1.3b) 


thus,  although  the  GO  fields  are  expressed  only  in  terms  of  the  axial  coordinate 
z.  the  transverse  coordinates  x  and  y  enter  implicitly  into  the  picture  through 


The  integrals  in  the  transport  equation  (1.2)  are  not  generally  amenable 

to  evaluation  in  closed  form.  This  is  simplv  because  the  z-dependence  of 

7*e  ,(z)  is  not  generally  known. 

a-1 


OBSERVATION 
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figure  1.  High-frequency  scattering  of  an  elec cronagnet ic 
field  from  a  smooth  convex  conducting  surface. 
The  dotted  lines  emanating  from  the  neighborhood 
of  the  specular  point  A  signifv  higher  order  ^ 

(k  '■)  contributions  to  the  observation  point  P. 


The  purpose  of  this  paper  is  to  overcome  this  difficulty.  Simple 
closed-form  expressions  for  a  generalized  GO  field  3t  the  observation 
point  ?  are  derived  by  combining  the  concepts  of  geometrical  optics  with 
the  diffraction  theory.  To  this  end,  the  expression  for  the  field 
amplitude  is  expanded  as  a  Tavlor  series  in  terms  of  the  coordinates  trans¬ 
verse  to  the  rav  direction.  It  turns  out  that  the  coefficients  of  this 
Tavlor  series  are  functions  of  the  radii  of  curvature  R,  and  R  at  A.  The 
field  amplitude  predicted  by  this  method  at  the  observation  point  P  is 
equal  to  that  oredlcted  bv  GO  augmented  by  a  correction  term  that  is  of 
order  '<  '  relative  to  the  GO  term. 

The  improvement  produced  bv  the  present  method  over  the  classical  GO 
method  can  be  explained  as  follovs.  The  GO  method  expresses  E1" (?) ,  the 
retlected  field  at  ?,  in  terms  of  that  at  a  single  point,  viz.,  the  specular 
point  A,  of  the  reflecting  surface  On  the  other  hand,  the  tedious  but  more 
exact  PO  Integral  expresses  £r  i  PI  in  terms  of  the  current  on  the  whole  surface 
of  _.  The  present  method,  however,  gives  pr (*)  in  terms  of  a  patch  of  current 
aro'ind  the  specular  point  A,  thereby  improving  the  accuracy  of  GO  without  paying 
the  high  nurerical  cost  of  rO. 

A  comparison  between  the  different  aoproximate  methods  used  in 
hlgn-f requencv  electromagnetic  scattering  is  given  in  Table  I. 

An  important  feature  of  the  present  method  is  that  it  can  be  associated 
with  the  launching  method  recently  developed  by  Mittra  and  Rushdl  [d]  that 
bvpasses  the  search  for  the  specular  point  A,  and  is  suited,  in  particular, 
to  situations  where  the  reflecting  surface  is  specified  onlv  numerically. 
Consequently ,  the  reflected  field  computation  with  generalized  GO  can  be 
carried  out  verv  efficiently. 


.OBjroDH  jaqjjjq 
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ON  THE  AMPLITUDE  VARIATION  OF  THE  00  FIELD 

This  section  presents  a  brief  summary  of  basic  GO  formulas.  It  also 
serves  as  a  clue  to  the  amplitude  dependence  on  transverse  coordinates  trvat 
will  be  adopted  in  Section  5. 

Tor  the  time  contention  exo  (j-^t],  a  fpica  1  scalar  component  of  a 
ray  optical  field  is  given  by  [5] 

f  "  A00(20)  DFfr.t^)  exp [- jks ( r) ]  (2.1) 

where  r  »  (x.v.i)  is  the  position  vector  in  rectangular  pencil  coordinates 
*  see  Fig.  2).  2  coincides  with  the  axial  rav  which  passes  through  the  origin 

and  is  measured  positively  ir.  the  direction  of  wave  propagation  from 

0  -  (0.  0.  z  ). 

s(r)  is  the  phase  function  given  bv 

s(r)  -  s (0)  ♦  2  +  l/2.-TQ  c  ♦  0 (  c  ')  (2.2) 


-  "  [x  y]‘  gives  the  transverse  coordinates  x  and  y  of  the  pencil- 
The  2x2  matrix  is  the  curvature  matrix  of  the  wavefront,  given  relative 
to  its  principal  directions  hv  the  diagonal  form  T„. 


0 


_ 1_ 


(2.3) 


AXIAL  RAY 


?i aire 


Asia,  and  paraxla.  ravs  in  a  douM”  Ji"«rgent  nenci 
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where  R^  and  R,  are  the  principal  radii  of  curvature  for  the  wavefront. 
DF(z,z.)  is  the  divergence  factor  of  the  pencil  given  by 


DFU.Zj) 


Pet  Q(z) 
Det  Q(zq) 


1/2 


1/2 


(W  (W 


(R^z) 


(R  ,-*-z ) 


(2.4) 


"or  convenience,  ve  express  the  divergence  factor  as  a  product  of  two  factors 


where 


OF  -  a0d0 


*0  " 


Vfo 

V* 


i /: 


(2.5) 


( 2. ba ) 


and 


(2.6b) 


a.  and  d^  an  l,e  viewed  as  one-*ii»ens lonal  divergence  factors  associated  with 

0  0 

the  principal  directions. 

Application  to  electromagnetic  waves  requires,  ^f  course, the  inclusion 
of  the  polarization  and, therefore. an  extension  of  (2.1)  to  a  vector  form.  The 
electromagnetic  fields  E(r)  and  H(r)  can  be  written  as 


E(r)  •  e(r )  exp  (-Jks(r) I 

Hlr)  •  h(r)  exp  (-jks(r) ] 

where  the  phase  function  Hr)  is  again  given  bv  (2.2)  and  the  amplitude 
vectors  e(r)  and  h(r)  are  related  via  Maxwell's  equations  bv 


(2.6) 

(2.7) 


V(h 


(2. 


h 


4 


J)1 


(2. 


where  z.  is  the  free-space  wave  impedance.  The  second  terms  in  (2.8)  and 
(2.9)  are  neglected  under  the  classical  GO  aporoxlnat ion. 

To  resolve  the  amplitude  vectors  e(r)  and  h(r)  into  longitudinal 
and  transverse  components  with  respect  to  the  paraxial  ray  through  r,  we 
need  to  introduce  three  local  orthonoraal  base  vectors  at  the  point  r.  Vi 
the  aid  of  (2.2)  and  (2. 3), we  find  the  unit  vector  in  the  direction  of  the 


paraxial  rav  as 


10 


Since  is  a  perturbation  ot'  the  axial  unit  vector  z,  its  component  ir.  the 

z  direction  is  given  to  the  order  of  o 1  \  while  its  components  in  the  directions 

_  2 

transverse  to  z  are  given  to  the  order  of  |o  onlv.  This  means  that  o  is 

3 

equal  to  unity  correct  to  the  order  of  c  3. 

Sow  we  construct  the  two  orthogonal  directions  and  transverse  to 
7s 


[1  -  1/2 


(R^x)- 


1/2 


I.R^  +  2)  (R,>2) 


z 


(2.11) 


•  -  1/2 


XV 


(R,^t)  (R 


777  x  *  (1  -  1/2  Tr^TT  1  y 


R+2 


(2.12) 


The  trlade  (?^,  J^)  satisfies  the  orthonorr.alit”  conditions: 


1/2 


(?k  ‘  V 


•  5u.  *  0(x;.) 

<■  1J 


k. 


1.  2.  3 


(2.13) 


where  •  is  the  icronecker  delta  and  0(x  )  denotes  terms  of  the  order  of 

*" v  with  -  ♦  *3.  For  the  special  case  when  r  is  on  the  axial  ray 

(x-0,  y-0).  the  local  base  vectors  (b  ,  -  ,  7s)  coincide  with  the  pencil 
base  vectors  (x,  y,  2).  At  anv  point  in  the  pencils  we  represent  the  amplitude 


vector  as 
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e(f)  -  Vl  «■  e232  ♦  «3a3 


(2.14) 


where  (e^,  e?)  are  the  transverse  components,  and  e.  Is  the  longitudinal 
component  of  e  with  respect  to  the  direction  of  wave  propagation.  Substituting 
(2.10),  (2.11)  and  (2.12)  into  (2.14),  we  have 


•  •*[<!-  1/2  «,  -  1/2  ** 


(R^z)*-'  1  *'  *  (R  +z)  (R,*z)  '2 


(Rj+t)  3 


*,  I 


♦  v  (  -  1/2 


XV 


(R ,*z)  (R,+z)  '1 


*  (1  -  1/2  (R^)2)  «2 


R  ,4-z  *3  ' 


#  r  x  y 

1  "  (R^z)  1  "  (R,*z)  *2 


(1  -  1/2  Tr^t:  -  1/2  >  *3] 


x  •  ’  ♦  y  el,  ♦  z  e^ 


(2.15) 


While  the  amplitude  of  the  scalar  rav  optical  field  (2.1)  is  independent  of 
the  transverse  coordinates  x  and  v.  Equation  (2.15)  reveals  that  the 


12 


amplitude  vector  *  is  a  function  of  the  transverse  coordinates.  If  we  ignore 
the  transverse  dependence  of  the  local  components  e^,  e,  and  e^,  then  the 
pencil  components  ej  ,  e',  and  e^  are  second  degree  polynomials  of 

— - —  and  — .  Using  this  result,  we  will  assume  a  generalized  ray  optical 

R-,+r 

field  of  the  fora  ; 

nrttv_2 

f„  -  (  *  A  (t)  xmyn)  exp  [-jk.s(r)l  (2.16) 

0  *■*  ®n 

m,n  •  0,0 

which  reduces  to  the  classical  ray  optical  field  when  x  -  y  -  0. 
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3.  FIELD  TRANSFORMATION  BETWEEN  TWO  PU.NAR  APERTURES. 

This  section  is  devoted  to  a  brief  presentation  of  some  notations 
and  definitions  which  are  needed  for  the  following  two  sections.  The 
reader  can  find  a  detailed  description  of  the  material  of  this  section  in  [6] 
and  [7|.  Let  f^p^)  represent  a  field  distribution  in  the  z  ■  plane,  whose 
transverse  coordinates  x^,  are  given  by  the  vector 


he  Fourier  transform  pair  is 


-T- 

F(t)  -  Ff  f (T) J  -  7  f(c)  e  '  do 


f(o)  -  f  (r(?)l  -  r(?)  e~JC  d? 


where  *  Is  a  two-dimensional  vector  in  the  spectral  domain,  given  by 


"CO  - 


and  the  superscript  T  on  a  vector  indicates  that  it  is  transposed. 


The  field  distribution  f,(c^)  in  the  plane  z  -  zy  due  to  an  aperture 
distribution  f Ac  in  the  plane  z  •  z  ^  (see  Fig.  3)  is  given  by  the  Fresnel- 
Xirchoff  integral  (6,  Eq.  3.5). 


fi(“i)  *  v(zI-zo)  ’ ‘  fo l- 0 ’  dc< 


where  •  is  the  free-space  wavelength,  k  "  —  is  the  free-space  wavenumber, 
and  R  is  the  distance  between  the  points  (.  2^)  and  („-  ,  s^),  given  by 


*"  '  <*i  -  V  *  'Vo  ‘  • 


Figure  1.  Coordinate  system  pertinent  to  the 

continuation  ot  a  field  distribution 
from  one  plane  to  another. 


for  all  Jq  and  of  Interest,  we  can  express  (3.4)  as 


R  -  U.-z,,)  ♦  - - r" 

1  0  2(«l-0) 


“VV 


“  T 

"0  “1 
<W  ' 


(3.5) 


(3.6) 


Substituting  (3.6)  into  (3.3)  and  waking  use  of  (3.1),  we  finally  arrive  at 


-Jk(«i-«0) 
f.(0.)  *  ‘'“77 - r - 

1  1  *<vV 

where  the  spectral  domain  vector  *  used  In  the  definition  of  the  Fourier 
transform  Is  given  by 


■j  2(*,-«n) 


f  <W 


e  2(xj-*o)  ) 


(3.7) 


(3.8) 


Equation  (3.6)  can  be  expressed  In  an  equivalent  form  of  a  Fresnel  transform 
[6],  but  the  above  form  is  sufficient  for  our  present  purposes.  In  the  next 


two  sections,  we  will  make  use  of  (3.7)  to  derive  the  GGO  formulas. 
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4.  APPLICATION  of  field  transformation  to  a  go  field. 

In  this  section,  we  will  show  that  if  a  ray  optical  field  exists 

at  the  aperture  plane  z  ■  z  ,  then  subject  to  the  Fresnel  approximation,  the 

o 

field  at  the  aperture  plane  z  ■  z^  will  be  that  ray  field  predicted  by 
geometrical  optics.  This  result  is  Important  in  its  own  right  and  serves  as 
a  preliminary  to  the  results  to  be  obtained  in  the  following  section.  From 
(2.1)  and  (2.2),  a  tvplcal  scalar  component  of  a  ray  optical  field  at  an 
aperture  plane  z  ■  z  ^  can  be  written  as 


vV  •  w  ex?  Mk  ^oHFoVo1 


(4.D 


m 

where  is  the  curvature  matrix  of  the  wavefront  at  the  plane  z  ■  z  . 
Substituting  (-.1)  into  (3.7)  vields 


W  • 


A, 


00(V 


-J  k  “1 


v(vV 


2(zrV 


-l/2oJ  ’T  To 


(4.2) 


where 


H  •  Jk(Q0  ♦ 


<VV 


I) 


(4.3) 


where  I  is  the  unit  matrix,  and  ^  is  a  complex  covariance  matrix  for  the 
zero-mean  two-dimensional  Gaussian  distribution 


2-  d.eirL  1/2 


The  Fourier  transform  of  the  above 


distribution  is  given  by  [8] 


Subst lcucing  (4.4)  into  (4.2)  and  making  use  of  (4.3),  we  finally  arrive 
at 


.  k  ,  -1 


(zj  e  1  e  2  1  *1  *0  °  1  0  1 


w  - 


00  0 


[det(Q0(r1-i())>I)) 


l /; 


(4.5) 


where  we  have  made  use  of 


(dec  H  l) 


-if: 


1/2 


i/: 


■  (dec  k  ) 


MW 


[  dec (Q 


(W  ♦ 


D) 


(4.6) 


We  now  proceed  to  show  c  it  expression  (4.5)  for  reduces  to  that 

predicted  bv  geometrical  optics.  In  view  of  (2.3)  ve  obtain: 


vw 


s. 


z 


z 


1 

0 


0 


0 


R,  ♦ 


(4.7) 


Subsequent lv ,  ve  get: 


[dec  (Q0(z1-2q)  ♦  I)] 


1/2 


'W  (W 
CW  (W 


(DF(zltz0)] 


(4.8) 


((Qrt  <*,-*«>  +  i>  -  I) 


(^-Zq)  ‘"*0  '‘1  ‘0 


S  ♦z 
1  1 


2  1 


'4.  9) 


where  is  the  curvature  matrix  of  the  wavefront  at  the  plane  z  •  z,  . 
Vith  the  aid  of  (4.3)  and  (4.1),  (4.5)  is  reduced  to 


1  -T  - 


’  Aoo(Io)  DF(Irx0)  exp  ,"Jk(*i4‘25i  Qri)1 


(-.10) 


which  is  the  field  distribution  predicted  by  GO.  It  is  interesting  to  note 
that  (4.9)  can  be  reduced  to  Equation  (3.23)  of  (5).  viz.. 


1  •— *  » 

qi  ‘  Jo  *  (W  1 


(4.11) 


which  is  the  standard  equation  for  phase  continuation  in  geometrical  optics. 
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5.  EXTENSION  TO  GENERALIZED  GEOMETRICAL  OPTICS  FIFLD. 

In  the  representation  (9.1)  of  the  ray  optical  field,  the  amplitude 
factor  aqq(2,))  is  3  function  of  the  axial  coordinate  z^  only.  In  this  section, 
a  generalization  of  that  field  will  be  considered,  where  the  amplitude  factor 
is  allowed  to  be  a  slowlv  varying  function  of  the  transverse  coordinates  and 
y().  In  addition  to  its  z-  dependence.  This  generalized  GO  field  which  was 
suggested  by  Eq.  (2.16)  is  repeated  here  for  convenience. 


w 


(  l 

®, n«0 ,0 


Aan(*0)  V”’  *Xpl*Jk(20"i  ?0  Vo)!  ' 


We  assume  that  the  field  (5.1)  exists  at  the  plane  ?^(z»z^)  which  is  normal 
to  the  reflected  ray  at  the  specular  oolnt  a  in  the  electromagnetic  reflecticn 

problem  of  Section  1.  (See  Fig.  9).  The  coefficients  A  (z  ),  which  are 

ir.n  0 

posslblv  complex,  can  be  expressed  in  terms  of  the  Incident  field  and  the  local 
surface  properties  at  the  point  A.  Our  task  here  is  to  use  the  results  of  Sections 
J  and  ■»  to  find  the  field  distribution  f^(c^)  at  the  plane  (z-z.)  normal 
to  the  reflected  r»v  at  the  observation  point  P. 

Substituting  (5...)  in  (3.6)  vields: 


(5.2) 
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where  Is  given  bv  Eq.  (-..3).  Equation  1.5.2)  expresses  fj(0,)  In  terms  of 

six  Fourier  transforms.  The  first  of  these  Is  given  bv  (4.4)  as 


«  »0 


-i/2?t  “  -i  7 


F(t)  -  F  {« 


■  1/2 
(jet  JJ  ) 


the  other  five  transforms,  obtained  in  Appendix  A,  are 


12  0  T  a  o0J  -|Vx,  *  1/2  lb»c)  V,  r(T )] 


(5.4«) 


F  vo • 


*  30  **  So]  *  fl/2  (b*c>  X.  *  d2v.J 


(5.ib) 


'  [x°  '  "  0  *  °]  '  [fir  (*i-  *o>  *  “S  ‘  1 


2  1  M"C  )  v,  )  r,-) 


(5.4c) 


r  (jo  e‘  ‘  ‘  0  "  :°]  "£fk(xi  “  V  *  (1  2  (b*c)  xi  *  dS)2] 


(5.4d) 


F[v°  V 


‘  a‘x.  *  1  2  (b*c)  v  )  (1  ;  (b*c) x. 


,  *  «%>] 


(5.  -«) 


where  the  non  dimensional  scalars  a“,  b* ,  c  and  d  are  given  bv 


%  (,rV  *  l) 


•  I  c  d2 

L  J 


rherer  ore,  equal  ion  (5.2)  for  the  field  distribution  at  tne  plane  ;  can 


rewritten  as 


f^U,)  •  DFii1>t0J  exp  [-]k  ( * r  ♦  1/2  p‘  ^  0^)1 


lzl"z0*  2  2 

A00  *  2j'k  ,;j  *20  *  Zd  *02  *  lb*CI  All  ’ 


♦  U*  A,  .  ♦  1/2  (b*<>  A  ,)  x 
1  1  011 


(1/2  (b+c)  A,  ♦  d*  A  ,)  v 

kJ  01  1 


4 

(a  A ,  ,  ♦  1  / »  (  b-*-c  )  A  ^ ,  ♦  1  /  2  a  ( b-H: )  A ^  ^  )  x ^ 


(1/4  (b*c)  A2q  ♦  d  A.,,  ♦  1/2  d‘  <b*c)  Aj  )  v‘ 


(1/2  a*  (b+c)  A,„  *  1/2  d‘  (b*c)  A„ ,  *  («‘d;  *  1/4  (b-*c)2)  A.,) 


Equation  (5.6)  is  the  GCO  formula  we  had  been  seeking  and  lc  has 
a  form  suitable  for  computer ional  work.  However,  additional  insight  into 
the  physical  interpretation  of  this  formula  can  be  gained  bv  writing  (5.6) 


in  a  special  form  as  follows.  If  we  choose  (Xg,  v^)  as  the  principal 
directions  (x.^,  Vgp) .  b  and  c  become  zero  while  a  and  d  reduce  to  a^  and 
d^  given  bv  Eqs.  (2.6a)  and  (2.6b),  respectively.  We  then  have 


l(riD>  *  VO  eXP  [-jk(Zl  *■  1/2  :1TD  5  id  ^  ID  * 


f a+n  <  2 

\  A  Am.n(V  (^X1D)  (d0  V 


a,n»0,0 


Z1_I0)  2D  2D 

- t -  <*q  A;o(V  *  d0  A02(Z0) 


We  show  in  Appendix  B  that  (5.6)  and  (5.7)  are  equivalent,  i.e.,  the  GGO 
formula  is  invariant  to  a  rotation  of  the  transverse  axes. 

Some  comments  on  the  GCO  formula  are  in  order: 

1)  Due  to  the  transitivity  of  the  Fresnel  transform  operator, 

the  expression  for  f^(c^)  in  terms  of  fg(Cg)  Is  self-consistent 
If  we  Fresnel-c rans form  either  fg(Cg)  of  Eq.  (5.1)  or  f^(e^)  of 
Eq.  (5.6)  to  obtain  f,(T,)  at  2  •  z,,  we  arrive  at  the  same 
result.  This  is  proven  in  Appendix  C. 

2)  The  classical  GO  amplitude  term,  viz.  A^g,  is  dominant.  For 

a  doublv  divergent  pencil  (R, >0,  R,»0),  this  term  decreases 


24 

with  the  pencil  advancement  as  does  the  divergence  factor 
a^dy  The  decav  of  the  higher-order  amplitude  terms  Is 
proportional  to  a^  d^  •  hence  this  decay  Is  higher  for 

a  larger  m. 

3)  The  amplitude  at  an  axial  point  Is  given  by  ( 5 . 8a > 


J 1  *  i  )  •>  n  ?  n 

yv— ^  (-oa2o(v  +  <wv> 


2  ^ 

Substituting  for  a*  and  d*  we  can  write  it  explicitly  as 


Aoo(to) 


1(W 


R,  -*-2  . 

[  (tTT~)  A°„(zJ  ♦  ( 


R  +z 


R  *■: 
*1  1 


20  0 


A02(Z0) ' 


(5.8b) 


The  0G0  amplitude  result  differs  from  the  GO  result  Ann(z0) 
bv  a  tern  of  order  k  For  large  and  finite  R1  and  R,, 
the  GGO  amplitude  expression  asvmptotically  approaches 


'  )k‘‘  "Ri*V  A2o(V  *  <W  am<V' 


00 


(5.9) 


4)  The  correction  term  in  (5.8)  can  be  obtained  from  the  transport  equation 
(1.2b)  for  m  -  1  if  we  have  the  following  two  conditions: 


a) 


>  > 


(5. 10) 
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b)  The  transverse  second  partial  derivatives  of  cQ  have 


the  following  z-dependence. 


5/: 


i  e 


,XD 


*D 


-  2  A  (z)  -  2  ( 

20  R+z 


R  ♦Z- 
- — 


1/ 

(— — 


7°  t 
A20  Z0 


■> . 

) ~e , 


>v. 


1/2 


i /: 


-o  ,  .  ,  ,R 1**0. 

A02  *  "  ‘  R ,  -*-z 


( ^ 


*0  ,  X 

A02  Z0 


The  tvpe  of  z-dependence  Indicated  by  (5.14)  is  the  same  as 
those  indicated  by  the  paraxlal-rav  analysis  (Eq.  (2.5))  and 
bv  the  Fresnel  transform  analysis  (Eq.  (5.7)).  As  pointed 
out  in  Section  1,  it  was  the  non-availability  of  the  z-dependence 
of  he  that  prevented  a  closed-form  evaluation  of  higher-order 


(5. 11a) 


(5.11b) 


GO  fields  at  the  observation  point. 
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APPENDIX  A 

COMPUTATION  OF  THE  FOURIER  TRANSFORMS  IN  EQUATION  (5.2) 


The  first  of  the  Fourier  Transforms  in  Eq.  (5.2)  is  given  by 
F(?)  -  F  (exp  (-1/2  H  ;0>  ) 


2s 


exp  ( 


-*T  •— I  -* 

-1/2  Tl  u  T 


T) 


(d«ti!  >1/2 

Using  (4.  3)  and  (5.5),  we  can  exoress  the  2«  2  matrix  ® 


—  1  V*0 

“  ’  }k 


a  b  ] 

* 

c  d“| 

>  J2 

L  J 

- 

so  that  (A.l)  can  be  rewritten  as 

„  ■>„  2  2  . 
F(*)  •  - - - ; — t  exp  (-1/2  (1  r.  ♦  (S+hrKn 

(det  H  > 

Sow,  rewriting  (3.1)  as 

. iff  - 

F(')  -  F(f(cQ)]  •  Jj  f(Cg)  e  do 

r  equivalent lv 

F(t.n)  -  f(x0,y0)  exp  (J(xqC  ♦  ygh))  dxQ  dvQ  . 


(A.l) 


(A.  2) 


(A.  3) 


(A.  4) 


we  can  arrive  at  expressions  for  the  last  five  transforms  in  Eq.  (5.2)  by 

2  2 

differentiating  (A. 3)  or  equivalently  (A. 4)  with  respect  to  t.  ' ,  i  ,  "  and 


respectively. 
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F  [xn  exp  (-1/2  oQ  H  oQ) 1  -  -J 


-  j/2  (2  C  +  (3+v)  'll  F(t) 


F  (y0  exp  (-1/2  cj  5  e0>]  -  -J  jF^,n> 


an 


J/2  [(S+>)  :  +  2a2  n)  F(t) 


F  (xj  exp  (-1/2  U  Oq> )  “ 


>~F(- .-) 
3  C2 


-  [a  -  1/4  (2aL  C  ♦  (2*>)  n)  )F(t) 


F  [y2  exp  (-1/2  cj  u  jq) )  - 


illili.il 

3n2 


-  (a  -  1/4  (  (3+-,)  5  ♦  2  a  n  )  ]  F(') 


F  l  Vo  exp  ('1/2  "o  ”  °o)]  - 


.  a*F(S.n) 

U5n 


-  [1/2  (S-^v)  -  1/4  (2a1C  ♦  (S+a)n)((3+Y)  ♦  2a2'')  ]  F(t) 


(A.  5) 


(A. 6) 


(A. 7) 


(A. 8) 


(A. 9) 


1 


I 
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(A. 12) 
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and 

—  L  — 1T 

la  ♦  a  1 


Zs 


b-t-c 


V. 


Equation  (A. 12)  is  equivalent  to  Cqs.  (A. 5)  and  (A. 6),  while  Eq.  (A. 13) 
is  equivalent  to  Eqs.  (A. 7),  (A. 8)  and  (A. 9). 


i 


♦  'i 
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APPENDIX  B 

INVARIANCE  OF  THE  GGO  FORMULA  TO  THE  ROTATION  OF  THE  TRANSVERSE  AXES 


Consider  the  case  when  the  transverse  directions  (x,y)  do  not  coincide 
with  the  principal  directions  (x^,  v^)  of  the  wavefront  but  make  an  angle 
with  respect  to  them  (Figure  5).  The  coordinate  vectors 


and  cD  - 

V  V. 


ire  related  bv 


-  T  p. 


-T  - 
■  T  o 


where  the  transformation  matrix  T  is  a  unitary  matrix  given  by 


S  C 


where  C  stands  for  cos?1  and  S  stands  for  sin?. 

Since  the  phase  quadratic  form  1/2  .  is  invariant  to  rotation 

of  the  transverse  axes  [5  ],  we  have 


1,2  h  Vo  '  1,2  ^  f  3*  ?T  « 


1/2  cT  Q, 


front  which  we  obtain 


•  9  •  «T 

Q0  "  ^  °0D  T 


(B.4) 


Using  (2.4),  (2.6)  and  (B.5)  we  rewrite  (6.21)  to  (6.23)  as 


b+c  •  -2  SC  V 


where 


(B.6) 

<  B .  7 ) 

(B.S) 


R;>r01  *0  d0 


3  <4 


so  that  (B. 6)  to  (B.S)  can  be  rewritten  as 


■>  i  ■)  '>  j 

d"  •  C*  ♦  S*  d~ 

o  o 


i  ■»  •>  2  1 

d*  •  S*  a*  -  C‘  d* 

o  o 


b  ♦  c  *  -SC  (a'  -  d“) 

o  o 


Next  we  relate  the  amplitude  coefficients  In  the  two  coordinate  systems. 
’Jsing  (5.1),  a  scalar  component  of  the  amplitude  vector  can  be  expressed  as 


-T 


ho  "  Aoo  *  Aio  A01  ^0  '  "0 


,D  .,D  ,D  ,  *0  *0 

A00  A10  A01  “0  *  ‘0 


'  20 

1/2  A 

4 

1/2  A 


11 


D 

11 


1/1  *u"! 

Aoi  J 


1/2  A 


4 


u 


*o 

°o 


IB. 10a) 


(B. 10b) 


Hence,  by  making  use  of  (B.D,  we  arrive  at 


A00  *  A00 


(B. 11a) 


A10  "  C  A10  *  s  A01 


(B. lib) 


A01  ‘  S  A10  *  C  A01 


(B . 11c) 


Au  -  2  SC  (A0,  -  A:f))  *  (C  -S‘)  A 


11 


(B.lld) 


A20  *  C*  A;o  *  S*  A02  *  cs  A 


11 


<  B. lie) 
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A 


:o 


2 

O' 


'02 


CS  A 


11 


CB. Ilf > 


Our  task  now  Is  to  show  that  the  scalar  amplitude  hj^  is  the  same  whan  expressed 
in  either  coordinate  svstem.  In  view  of  (5.6)  and  (5.7),  we  want  to  show  that 


h°  -  DF 


f  *+nc_2 
\n,n“0,0 


Ann'V  (VlD)  (d0V10> 


U.-t  )  2  -j 

-  >  “o  A20U0>  *  j 


o  4<V  } 


(8.12) 


and 


( 


(VV 


h  •  DF  {  A  -)  - r - (  *‘A20  ♦  d"A02  ♦  1/2  (b*c)  An  > 


♦  {•  Aio  *  1/2  (b4°  A01‘  X1 


[1/2  (b*c)  Al0*  d‘A01]  yx 


[a"  A,0  *  1/*  (b+c)‘  Aq;  ♦  1/2  a‘  <  b+c )  A^]  x‘ 


(1/e  (b^c ) “  A  *  d  A  +  1/2  d  <b*c)  Au]  y‘ 


(1/2  a"  (b+c)  A ,g  *  1/2  d‘  (tr*-c)  A ^ ,  ♦  (e‘d‘  *  1/4  (b+c)“)  A^] 


v} 


(B.1J) 


3b 


are  equal. 

Substituting  (B.2)  and  (B.ll)  in  (B.12)  we  get 


ft  C  <VV  7  1  2  1 

»;  •  DF  (_*oo  -1  -V-  1  ‘*5  c'  *  Jo  s‘>  A:o 


y  y  y  2 


tJ0  s  *  40  L  1  *02  +  Cs  ‘*0  "  V  *ll' 


*0  (C  xl  ♦  s  V  (c  *10  *  s  *01 } 


d0  (-s  xi  *  c  V1>  (-s  *10  *  c  *01* 


y  y 


*  *g  <c  *1  *  s  y^*  ^c‘  A,0  •*■  S*  Aq2  f  CS  Au> 


♦  J0  (-S  xl  *  C  V  (s*  *20  *  C‘  *02  "  CS  *11^ 


♦  «„  dj  <c  »!  *  s  «,)  C-S  X,  *  C  V 


jiu  SC  (*0,-»,0>  *  <C:-S:)  (»•!' 


which  reduces  to  <B.13)  bv  virtue  of  (B.b')  to  (B.8').  This  shows  that 
representing  the  GCO  formula  relative  to  the  principal  axes  (x^,  y^)  is 
equivalent  to  representing  it  relative  to  anv  other  set  of  orthogonal  axes 


(x.  v). 
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AP  PEND IX  C 

TRANSITIVITY  OF  THE  CGO  FORMULA 


Due  to  the  transitivity  of  the  Fresnel  transform  operator  [6],  the 

expression  for  f^(o^)  In  terms  of  f^fp^)  Is  self-consistent  In  the  sense 

that  If  we  Fresnel-transf orm  either  f^(pg)  ot  Eq.  (5.1)  or  fj(c^) 

-♦ 

Eq.  (5.7)  to  obtain  f.(c,)  at  *"*■»•  ue  arrive  at  the  same  result.  In 
view  of  the  results  In  Appendix  B,  we  lose  no  generality  In  the  following 
proof  If  we  take  the  transverse  coordinates  as  the  principal  ones. 
Starting  from  f and  using  Eq .  (5.7)  we  can  write 


fl(V  •  *01  d0l  *X?  (-Jk  (zi  *  1/2  ®U>  51D  •rtD)1 


f 


x+ry-2  D  7  «  ? 


\  .  Amn  (I0)  <*0i  X1D)  (d0l  'V 

lx.n-0,0 


^ 1  a “  *  -)  ^  2  D  2  D 

k  *01  A20  (Z0)  *  d"i  A02  ‘‘o' 


(C.  1) 


1  •  1,2 


where 


v*0 

,  -  — a - 2 

01  V*i 


1  -  1.2 


and 


(C. 2a) 


s  I*-.**, 

1“  »  — - 

0i 


(C. 2b) 


1  -  1 
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Starting  from  f^(„'^)  ot  E8-  (C.l)  ,  ve  can  rewrite  a* 

1/2 


f,(o,)  -  ( 


(Rj+tj)  (R,«-2,) 


(R^*,)  (R,+2, 


J  1  *01  d01  *XP  '-Jk  “l*  1/2 


{■ 


nrt-ru.2 

T 


a,n«0,0 


D  .  ,  2m 

A  (2  )  «n 

mn  G  ’J 


*  <2n  f  V!l  x  1  fV!i  i 

1  0  1  [_R1+*;  “J  LK2+*2 


zr20  ,  2  .D 


j  k  l*oi  a2o  (V  *  d01  Ao:  (V  1 


e2~*  i  r  vii  ad 

k  LR^2  20 


**2*zl  D 

(s,)  ♦  ~ —  A.,  (2  ) 


T  R,fr,  02  1 


which  reduces  to 


-T  - 
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m-*-tv^2 
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^‘>4’*0  2  D 
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) 


♦  (  (2-2,)  *  (2,-2,) 
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To  show  the  agreement  of  Eq.  (C.4)  with  Eq.  (C.2)  for  1-2,  we  need  only 
to  consider  the  k  ^  term,  which  Is 


(VV  (Rl"V  *  (z2"ll)  (Rl>Z0) 

(R^z,)  (R^z^ 


<W  A20  (V 


<VV  (R:+22)  +  (*2~*l)  <V*0)  n 

(  ~  T*2+t'2)  'R^+Zj)  )  (R2+*0)  A02  (i0)  1 


j(z:-z0)  R,+z 


1‘0  ,D  .  .  R2*Z0  D 

Rj+z,  A20  *0  *  R,+z.,  A02  (V  ^ 


(W 


R  f  *02  A20  (V  *  d02  A02  (V  ! 


(C. 


The  nroof  of  transltlvltv  of  OCO  formulas  Is  complete. 


40 


REFERENCES 

[1)  W.V.T.  Rusch,  "Reflector  Antennas, "  in  Numerical  and  Asymptotic 
Techniques  in  Electromagnetics,  R.  Mittra,  Ed.  New  York: 

Spr  inger-V'er  la*.  1975,  pp.  217-256. 

[2)  S.  U.  Lee,  "Electromagnetic  Reflection  from  a  Conducting  Surface: 
Geometrical  Optics  Solution,"  IEEE  Trans.  Antennas  Propagat.  , 
vol.  AP-2J,  no.  2,  pp.  134-191,  March  1975. 

[3)  S.  W.  Lee,  "Private  Communication." 

[4)  R.  Mittra  and  A.  M.  Rushdi,  "An  Efficient  Approach  for  Computing  the 
Geometrical  Optics  Field  Reflected  from  a  Numerically  Specified  Surface, 
IEEE  Trans.  Antennas  and  Pro?.,  vol.  AP-27,  no.  6,  November  1979. 

[5)  S.  W.  Lee,  "Geometrical  theory  of  Diffraction  in  Electromagnetics, 

Vol.  1:  Geometrical  Optics,"  Electromagnetics  Laboratory  Report  No.  78 
Universitv  of  Illinois  at  Urbana-Champaign,  Mav  1978. 

[6)  R.  Mittra  and  P.  L.  Ransom,  "Imaging  with  Coherent  Fields,"  Proceedings 
of  the  Symposium  on  Modern  Optics.  Polytechnic  Institute  of  Brooklyn, 
March  1967,  pp.  619-647. 

[7)  W.  T.  Cathev,  Optical  Information  Processing  and  Holography.  New  York: 
John  Wiley  4  Sons,  1974. 

[81  A.  Papoulis,  Probability,  Random  Variables  and  Stochastic  Processes. 

New  York:  McGraw-Hill,  1965. 

.  D.  Luenberger,  Optimization  by  Vector  ,i2i£JL  Methods .  New  York: 

John  Wilev  4  Sons,  1969. 


